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Abstract 

The large mixing angle (LMA) MSW solution to the solar neutrino 
problem seems favored by the data at the moment over the small 
mixing angle (SMA) MSW solution and the vacuum (VAC) solution. 
In this paper the various main types of models of neutrino masses and 
mixings are studied from the point of view of how naturally they can 
give the LMA solution. Special attention is given to a very simple 
type of "lopsided" SU (5) model. 
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1 Introduction 



The main solutions to the solar neutrino problem are the SMA solution 
(small mixing angle MSW), the LMA solution (large mixing angle MSW), 
the LOW solution, and the VAC solution (vacuum oscillations). The exper- 
imental situation has been and remains very ambiguous. However, recent 
results seem somewhat to disfavor the SMA and VAC solutions. In fact the 
superKamiokande collaboration claims that they are excluded at the 95% 
confidence level 

On the other hand, a survey of the hundreds of models of neutrino masses 
and mixings published in the last three years shows that most of them yield 
the SMA or VAC solution, and even some that claim to obtain the LMA 
solution are only marginally consistent with the latest global analyses of the 
data. The purpose of this paper is to look at the main types of models of 
neutrino masses and mixing angles that have been proposed in the literature 
from the point of view of their ability to yield the LMA solution in a natural 
way. There are two aspects of this question that can be distinguished. First, 
one can ask whether a certain scheme or model can fit the LMA solar solution 
with some choice of model parameters that is not too badly fine tuned. Sec- 
ond, one can ask whether the model explains the LMA values of the neutrino 
masses and mixing angles. In order to say that a theoretical model really 
explains them, something close to the LMA best-fit values should emerge 
automatically when the parameters of the theoretical model take their most 
"natural values" . If a model parameter that is a priori of order one must 
be set to a value of ten or a tenth in order to fit the neutrino masses and 
mixings, one has accomodated them but not really explained them. What 
our survey will show is that of the great number of models that now exist 
in the literature, few can be said to provide an explanation (in this sense) of 
the LMA values of tan 2 9 so i and Am^. 

What are the values that are to be explained? A recent global analysis 
of Bahcall, Krastev and Smirnov 0] (published before the SNO results Pj) 
gives as the best fit values for the LMA solution Am^, = 4.2 x 1CT 5 eV 2 , 
and tan 2 # so ; ~ 0.26. However, the allowed region is fairly broad. The 90% 
confidence-level allowed region given in || extends in Am 2 o( from about 
2 x 10^ 5 eV 2 to 10~ 4 eV 2 , and in tan 2 9 so i from about 0.15 to 0.5; i.e. about 
a factor of two in either direction for both parameters. A significant aspect 
of the fits in p[ is that they nearly exclude exactly maximal mixing for the 
LMA solution. The 95% confidence level contour extends in tan 2 6 so i only up 
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to 0.55, and the 99% contour up to 0.7. A more recent analysis by Bahcall, 
Gonzalez-Garcia and Pena-Garay [|J, which includes the SNO results, gives 
similar results for the LMA solution; their best-fit value of tan 2 9 so i being 
about 0.4. However, the allowed LMA regions of |4j] extend somewhat higher 
in tan 2 9 so i than those obtained in || . For example, the fit of Fig. 9 of Q , 
in which sterile neutrinos are excluded a priori, gives a 95% confidence- level 
region that extends up to about tan 2 ^ so ; = 0.8, and a 99% confidence- level 
region that extends slightly past tan 2 6 so i = 1. 

As we shall see, the issue of how close to exact maximal mixing tan 2 8 so i 
is allowed to be is crucial for deciding whether several kinds of models can 
naturally give an acceptable LMA solution. In this paper we shall say that a 
model gives an LMA solution that is in comfortable agreement with the data 
if it predicts tan 2 6 so i < 0.8. It is convenient to express the mass-squared 
splitting Am 2 oi as a fraction, which we shall call r, of the mass-squared split- 
ting relevant to the atmospheric neutrino oscillations: r = Am 2 D ;/Am 2 im . 
The best fit values give r ~ 1.4 x 10~ 2 . 

The rest of the paper is organized as follows. In section 2, we shall look 
at the basic non-see-saw approaches to neutrino mass. We shall see that 
models based on such approaches have been constructed which can fit the 
LMA solution, but for the most part not comfortably, either because r tends 
to come out too small or because the solar mixing angle tends to come out 
too close to maximality. In other words, most of the non-see-saw models 
that can fit the LMA solution do not really explain it in the sense that we 
have defined. In section 3, we look at see-saw models. Here too, most of 
the published models give the SMA or VAC solutions. However, we show 
that there are some reasonably simple "textures" that can reproduce nicely 
the LMA values of the neutrino masses and mixings. However, it remains 
unclear whether these simple textures can arise in simple models. 

In section 4, we look at a well-known model that is particularly interesting 
for two reasons: (a) it is very simple in conception, and (b) it can explain at 
least the LMA value of the solar neutrino angle, although it does not explain 
the value of the neutrino mass splitting. It is an SU(5) grand unified model 
with an abelian family symmetry. We shall analyze this model in some detail 
both analytically and numerically. We shall show how the predictions of this 
model can be studied statistically in a completely analytic way by assuming 
that the unknown parameters of the model have Gaussian distributions. This 
method should be easily aplicable to many other kinds of models. 

Section 5 is a brief summary. 
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2 Non-see-saw models 



A. Non-see-saw models where 9 atm comes from M v . 

In non-see-saw models the mass matrix M v of the three light neutrinos is 
typically generated by new low-energy physics. It therefore has no relation, 
or only a very indirect relation, to the Dirac mass matrices of the charged 
leptons, the down quarks, and the up quarks (which matrices we denote L, 
D, and U respectively). This has the great advantage of making it easy 
to explain why the atmospheric neutrino mixing angle is very large (Uys, = 
sin6 atm = 0.7) while the corresponding quark mixing is so small (V c b — 
0.04). If the Dirac matrices are assumed to be hierarchical, then they would 
naturally give the small mixing angles seen in the quark sector. But if M v 
is unrelated to the Dirac mass matrices, it could easily have a very different 
form with large off-diagonal elements that gives large mixing angles. Non- 
see-saw models based on this idea are called Type 1(1) in p. 

The tricky question for this type of model is to explain why Am^ <C 
Am 2 atm . If the large mixing comes from diagonalizing the 2-3 block of 
M u , one would expect that m 2 and m 3 , the second and third eigenvalues 
of M v , would have similar magnitudes, in which case typically so would 
Am 2 sol = m 2 — m\ and Am^ tm = m\ — m\. The challenge then is to reconcile 
the hierarchy seen in the eigenvalues of M v with the large atmospheric mixing 
angle. To do this requires a special form of M v . Two special forms have 
been found viable in constructing realistic models, one leads to a so-called 
"inverted hierarchy" m\ = m 2 ^> m 3 , and the other to the ordinary hierarchy 
mi <C wi2 "C ?7i3 ■ We shall consider these in turn. 

Inverted hierarchy models. 

Inverted hierarchy models have the following special form for M v : 



Here c = cos# and s = sin#, where 9 ~ 1, and -C M. One can 
diagonalize this matrix in stages, the first step being to rotate by angle 9 in 
the 2-3 plane, bringing the matrix to the form 
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One sees immediately that m\ = m 2 = M ^> m 3 . The mass-squared splitting 
relevant to atmospheric oscillations is Am 2 tm = M 2 , whereas the splitting 
relevant to solar oscillations is Am^ = 2(mu + m' 22 )M, which is much 
smaller, as required by all the viable solar solutions. The atmospheric angle 
gets a contribution 6 ~ 1 from diagonalizing the 2-3 block of M v , so in 
the absence of some unlikely cancellation it will be large, as observed. But 
what of the solar angle? From the fact that the 1-2 block of Eq. (2) has 
a pseudo-Dirac form it is apparent that the solar mixing angle will be close 
to maximal. Consequently, inverted hierarchy models cannot give the SMA 
solar solution, but rather give "bimaximal" mixing. 

The inverted hierarchy form of Eq. (1) can arise in several plausible 
ways. One example is the Zee type of model ||. In the Zee model |7| 
there is a singly charged singlet scalar field h + , which is allowed by the 
standard model quantum numbers to couple (antisymmetrically) to both a 
pair of lepton doublets (h + LiLj) and a pair of Higgs doublets (h + Q> a $b), 
assuming that more than one Higgs doublet exists. If both types of coupling 
are present, a conserved lepton number cannot be consistently assigned to 
h + , and consequently AL = 2 Majorana masses for the left-handed neutrinos 
arise at one-loop level. The resulting one-loop mass matrix has the form 



For c <C a ~ 6, this has the desired inverted hierarchy form. 

The inverted hierarchy form can also arise in models with an approxi- 
mately conserved L e — — L T lepton number. If this quantum number is 
exactly conserved, then only the 12, 21, 13, and 31 elements of M v can be 
nonvanishing. If there are small violations of L e — — L T the form in Eq. 
(1) can result M. 

The question of present interest to us is whether the inverted hierarchy 
can give an acceptable LMA solution. To answer this one must look more 
closely at the solar neutrino mixing angle. This is given by 6 so i = Q\ 2 — 9f 2 , 
where the two angles on the right-hand side are the contributions that come 
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from diagonalizing M v and L respectively. From Eq. (2) it is easily found that 
tan 2^ 2 = 2M/(m 22 -m n ), so that 9\ 2 ^ 7r/4 - (m' 22 - m u )/4M. We have 
already seen that r = Arn% ol / Am% tm = 2{m' 22 + m\i) / M . Requiring that this 
be of order 10 -2 as needed for the LMA solution, and assuming that there are 
no accidental cancellations, one has that 9\ 2 = 7r/4+O(10~ 3 ). If 9\ 2 vanished, 
this would give tan 2 9 so i = 1 + O(10 -3 ), which is too close to maximal mixing 
to be in comfortable accord with the global fits. However, one expects that 



relative to the contribution from M v , and can therefore increase or decrease 
tan 2 8 so i from unity. If one assumes that 9\ 2 = 0.07 and has a relative minus 
sign to 9" 2 , then tan 2 9 so i = 0.75, which is consistent with the global LMA 
fits. However, one can see that the tendency of inverted hierarchy models is 
to give solar mixing that is closer to maximality than to the best-fit LMA 
value of tan 2 9 so i ~ 0.4. This is one reason why many of the published 
inverted hierarchy models claim a better fit to the VAC solution than to the 
LMA solution 0. A significant reduction of the experimental upper limits 
on tan 2 8 so i would make the inverted hierarchy idea much less plausible as 
an explanation of the LMA solution. For example, a value of tan 2 6 so i = 
0.5, would imply in the inverted hierarchy context that tan#f 2 — 047 = 



2.5wm e /m M , which would require a very special form of the 1-2 block of L. 



Ordinary hierarchy models. 

The other possibility for non-see-saw models that gives a large atmo- 
spheric neutrino mixing angle coming from M v and a hierarchy in the mass- 
squared splittings is 



Here, again, c = cos# and s = sin#, where 9 ~ 1, and m^- M. As written, 
the 2-3 block of the matrix has vanishing determinant; however, it is assumed 
that there are small corrections to these elements, which we have not written. 

As in the case of the inverted hierarchy models, one can diagonalize this 
in stages, starting with a rotation by angle 9 in the 2-3 plane. The result of 
such a rotation is to bring the matrix to the form 




This contribution can have any complex phase 
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(5) 
\ m' 13 M 

Because of the small corrections to the 2-3 block that were just mentioned, 
the 22 element in Eq. (5) does not vanish, but is small compared to M. This 
matrix gives Arn^j = O(m^) and Am^ m = M 2 . Thus the right hierarchy of 
splittings for any of the solutions can be achieved for the appropriate values 
of rriij/M. In contrast to the inverted hierarchy form, this form can give 
either small or large 6 so i, and in the large- angle case there is no preference 
for values of 8 so i that are very close to maximal. 

The form in Eq. (4) is clearly special in the sense that the 2-3 block is 
approximately of rank one. This would be unnatural unless some symmetry 
or mechanism guaranteed it. One possibility is that this form arises from a 



nonabelian flavor symmetry pOfl , however, this is difficult to achieve. Rather, 
almost all published models that achieve this form in a natural way use 
the idea of factorization. The idea of factorization is that the dominant 
contribution to the neutrino mass matrix has the form (M v )ij = fifj, which 
is obviously of rank one. If f\ <C f 2 , fs, this dominant term reproduces the 
large elements in Eq. (4). The condition that f\ is small compared to f 2 and 
fz is necessary to satisfy the experimental constraint that U e ^, < 0.15. One 
drawback of most models based on factorization is that they do not explain 
why fi is small. 

A factorized form can arise in various ways in non-see-saw models. A 
much studied example is supersymmetry with terms in the superpotential 
that violate both lepton number and R-parity. Cubic terms of this type are 
\ij k LiLjE£ and X^ k LiQjD\. The latter leads to one-loop AL = 2 neutrino 
mass diagrams, in which a neutrino converts into a virtual quark-squark pair. 
Assuming that the LR squark masses are proportional to the corresponding 
quark masses, this diagram gives (M u )ij oc \'* kl \'j lk md k rnd r Consequently, the 
b-quark/b-squark loop dominates, and gives a contribution that is propor- 
tional to A'jggA^gmf, which obviously has a factorized form. This gives only 
the heaviest neutrino mass, m 3 . The second largest neutrino mass comes 
from a similar diagram with both b and s quarks/squarks in the loop. Con- 
sequently, one has that r = (m 2 /m^) 2 ~ (m s /mb) 2 ~3x 10~ 4 . This is much 
smaller than the value of 1.4 x 10 -2 preferred by experiment; however, there 
are several unknown parameters that come into this calculation, such as the 
couplings X' i j k , so that nothing prevents the right LMA value of r from being 
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obtained [IT]. However, the model does not really explain the magnitude of 
r. 

We have only considered the effects of the cubic lepton-number-violating 
and R-parity-violating terms in the superpotential. There are also in general 
bilinear terms of the form LiH u . These have the effect of mixing leptons 
and Higgs fields, and so allow the sneutrino fields to acquire non-vanishing 
vacuum expectation values. That, in turn, through the sneutrino- neutrino- 
neutralino coupling gives a tree-level neutrino mass in which the neutralino 
plays the role of "right-handed neutrino". It is easily seen that this tree- 
level mass has a factorized form and gives mass only to one neutrino, i.e. 
777.3. The other neutrino masses, m 2 and m 1; arise from the one-loop di- 
agrams previously discussed. In consequence, in such models where both 
cubic and bilinear R-parity-violating terms contribute to M v one expects 
that r = (m2/m^) 2 ~ (loop/tree) 2 <C 1CT 2 . For this reason, most analy- 
ses of supersymmetric models in which the bilinear R-parity-violating terms 
contribute to M v conclude that there is much more parameter space for the 
VAC solution than for the LMA solution, i.e. the LMA solution requires 
special choices or tuning of parameters fl2fl . However, in [|l^] it is shown 
that under certain assumptions (specifically, that there are only bilinear R- 
parity-violating terms and that the SUSY-breaking terms are non-universal) 
the LMA solution can be achieved without fine tuning. Nevertheless, it 
seems, on the whole, that the SUSY models with R-parity breaking do not 
do well in explaining the LMA value of Am 2 ol . 

Another possibility for obtaining an approximately factorized form for M v 
that has been much studied in the literature is called "single right-handed 
neutrino dominance" (SRHND) |L4|]. As the name suggests, the idea here is 
that instead of there being three right-handed neutrinos, one in each family, 
as there are in typical grand unified theories or Pati-Salam models, there is 
just one right-handed neutrino, N c , which can have mass terms MrN c N c + 
fi(uiN c )(H). Integrating out N c gives a rank-1 factorized contribution to 
M v . If one assumes that N c couples with almost equal strength to the /x and 
r neutrinos, and (for some reason not generally explained) only weakly to 
the electron neutrino, the large terms in Eq. (4) are reproduced. 

One way to explain the smallness of the coupling of N c to the electron 
neutrino would be to impose a symmetry that distinguishes v e from the 
v yi and v T (but does not distinguish the latter from each other). Such a 
symmetry would also tend to suppress mixing between the v e and the heavier 
neutrinos, and thus give the SMA solar solution, as in the model of [15]. 
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Since the right-handed neutrino only gives mass to one neutrino, some 



other mechanism must be found to give mass to the other neutrinos. In [16 



the lighter two neutrino masses arise from loop effects. In []T7| they arise at 
tree level from integrating out other heavy states that have different flavor 
quantum numbers than N c . In [18| they arise from operators of the form 
UiUjH u H u /Mpi, which it is argued are generally there anyway, the idea being 
to avoid having to invent new beyond-the-standard-model physics to account 
for each kind of neutrino mass. In all these cases, mi and vri2 are much less 
than 7713, though the specific reason is different in each case: in [16| they are 
suppressed by loop factors, in [IT] by small flavor-breaking parameters, and 
in [HJ by Mp/Mpi. That SRHND models tend to give a strong hierarchy in 
neutrino masses is what one would naturally expect. Since the mechanisms 
that generate the largest neutrino mass and the other neutrino masses are 
different, there is no reason a priori that they should yield masses of similar 
scale. Rather, it would be a coincidence calling for an explanation if they did. 
Because most SRHND models give m 2 <C m 3 they yield the VAC solution or 
SMA solution to the solar neutrino problem rather than the LMA solution 



16, 17, 18 



To obtain the LMA solution, one wants 777.3 and to be only about a 
factor of ten in ratio. This suggests that they arise from the same basic mech- 
anism. One possibility is that all three neutrino masses arise from integrating 
out right-handed neutrinos, but that one of those right-handed neutrinos is 
somewhat lighter than the others and so dominates to some extent, but not 
by a large factor. However, this would really be just a special case of the 
ordinary see-saw mechanism, which we will discuss in the next section. In 
fact, the structure given in Eq. (13) is really based on this idea, which was 
proposed in [IHJ . 



B. Non-see-saw models where 6 atm comes from L. 

There is another class of non-see-saw models in which the large atmo- 
spheric neutrino angle comes predominantly from the diagonalization of the 
charged lepton mass matrix L. (This class of models is called Type 11(1) in 
||.) This has the advantage that it becomes easy to reconcile the largeness 
of atm with the smallness of Am^/Am^, since the former comes from L 
while the latter comes from M v . On the other hand another issue arises for 
this class of models, namely explaining why the CKM angles are small. Since 
the form of L is such as to give a large mixing angle 9 atm , one would naturally 
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expect that the Dirac mass matrices D and U of the quarks would be such 
as to give similarly large contributions to the CKM angles. The point is that 
in many kinds of models the Dirac mass matrices L, D, and U are closely 
related to each other. 

One possibility is that there are indeed large contributions to the CKM 
angles coming from U and D, but that these nearly cancel. This possibility 
is realized in a much-studied class of models based on the idea of "flavor 
democracy" [p0~ |. In flavor democracy models it is assumed that all the Dirac 



mass matrices have approximately the "democratic" form 

(6) 




This form can be enforced by permutation symmetries among the three 
families. In the limit of exact flavor democracy, the matrices U and D are 
exactly of the same form, so that flavor mixing in the quark sector cancels 
out. On the other hand, it is assumed that the neutrino mass matrix M v 
has a very different form. In most papers it is assumed to be approximately 
proportional to the identity matrix, though in some papers it is only assumed 
to be nearly diagonal. As a result, for the leptonic mixing angles there is no 
cancellation such as makes the CKM angles small. 

The MNS mixing matrix for the leptons is given by Umns — U\U V , where 
Ul and U v are the unitary matrices that diagonalize respectively L^L and 
MlM v . If L has exactly the democratic form, then 

/ l/y/2 -1/V2 \ 

ul = l/VE i/Ve -2/Ve = u FD . (7) 

V 1A/3 1/V3 l/y/3 I 

If the mass matrix of the neutrinos is exactly diagonal, then Umns — U\ = 
Ufd- This would give sin 2 29 atm = 8/9, which is consistent with the data, and 
tan 2 6 so i = 1, i.e. exactly maximal mixing for solar neutrinos. However, the 
matrix L clearly cannot have exactly the democratic form, as that is rank one 
and would give m e = = 0. There must therefore be small corrections to L 
coming from the breaking of the permutation symmetries. These corrections 
not only generate masses for the electron and muon but also make the angle 
6 so i deviate from maximality. For the simplest and most widely assumed 
form of these corrections to L, one can calculate the corrections to 8 so i in 
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terms of yjm e /m^. One finds, still assuming that M v is diagonal, that 

4 , 

tan 2 6 sol = 1 - -p^jmjm^ = 0.84, 

or equivalently sin 2 29 so i = 0.993. This is too close to unity to be in com- 
fortable agreement with the LMA global fits. Almost all published models 
based on flavor democracy have tan 2 6 so i = 1, or else obtain smaller values 
by fine-tuning. However, Tanimoto, Watari, and Yanagida have a version 
in which there are small corrections to M v that can reduce tan 2 6 so i to the 



region preferred by the LMA fits [21]. While this shows that it is possible 



within the flavor democracy framework to construct LMA models that can 
fit the data, it does not appear that flavor democracy does a good job of 
explaining the LMA value of tan 2 # so ;. Flavor democracy is more naturally 
compatible with the VAC or LOW solutions. 

We may summarize the situation by saying that most schemes that have 
been proposed based on non-see-saw mechanisms neither very comfortably fit 
nor really do much to explain the values of the neutrino parameters required 
for the LMA solution to the solar neutrino problem. The great majority of 
non-see-saw models in the literature more naturally give the SMA or VAC 
solutions. There are exceptions, which we have noted above. How close 
tan 2 9 so i is to 1 is a crucial issue. 



3 See-saw models 

The see-saw mechanism is usually associated with grand unification. In 
S'O(IO) grand unified models, and in most other unified schemes except 
SU (5), the existence of one right-handed neutrino for each family is required 
to make up complete multiplets of the unified group. Moreover, the see-saw 
formula M v = — iV^M^iV, where Mr is the Majorana mass matrix of the 
right-handed neutrinos and N is the Dirac mass matrix of the neutrinos, 
gives neutrino masses in the range required by experiment if the scale of Mr 
is near the grand unified scale. Thus both the existence of neutrino masses 
and their magnitude are elegantly accounted for by the related ideas of grand 
unification and the see-saw mechanism. In this section, we shall therefore 
assume that we are dealing with a grand unified model. 

A. See-saw models where 9 atm comes from M v . 
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In models based on 50(10), there is generally a close relationship among 
the four Dirac mass matrices N, U, D, and L. Indeed, in the minimal 50(10) 
model (which is too simple to be realistic) N = U oc D = L. The smallness 
of the CKM angles and the small interfamily mass ratios of the quarks can be 
explained by assuming that the matrices U and D are "hierarchical" in form. 
There are two simple kinds of hierarchical mass matrix that are frequently 
encountered in models 



where e' <C e <C 1- The entries in these matrices as written are to be 
understood as giving only the order in the small parameters of the entries. 
The first form in Eq. (8) has what may be called a "geometric hierarchy", 
since an off-diagonal element is of the same order as the geometric mean 
of the corresponding diagonal elements. The second form in Eq. (8) has 
what may be called a "cascade hierarchy", since the matrix is made up of 
successive tiers, a 1-by-l, a 2-by-2, and a 3-by-3, of ever smaller magnitude. 
Both forms in Eq. (8), if applied to the quark masses, give V& ~ e, V us ~ e'/e, 
and V ub ~ e' ~ V us V cb . 

While there is as a rule a close relation among the four Dirac mass matri- 
ces in 50(10), the Majorana mass matrix Mr of the right-handed neutrinos 
can be quite different in form. For example, in minimal 5O(10) the Dirac 
mass matrices all come from the same term, 16jl6jl0#, whereas the matrix 
M R comes from different terms, either 16jl6,,126# or 16jl6jl6#16#. A 
reasonable hypothesis is that the CKM angles are small because of the hier- 
archical nature of the Dirac matrices, while the largeness of 9 a tm and possibly 
of 9 so i has to do with the very different form of Mr. Models based on this 
idea were classified in || as Type 1(2). 

A potential difficulty with this idea is that if the Dirac mass matrix of the 
neutrinos N has a hierarchical form it tends, through the see-saw formula, to 
make M v also have a hierarchical form, indeed a more strongly hierarchical 
form. For example, suppose N = diag(e',e, 1), and we parametrize Mr as 
{M^ij = ciij. Then 






(9) 
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If all the djj are of the same order, then tan 28^ — 26023/(033 — e 2 «22) ~ £• in 
order for 8 atin to come primarily from diagonalizing M u , one must have 8^ ~ 
1. Clearly, this is only possible with some special form of Mr. One possibility 
is that 023/033 ~ e -1 . If this is true, then there is a hierarchy among the 
elements of Mr that is related to the hierarchy among the elements of N. 
That is, the atmospheric neutrino mixing angle is of order unity because 
of a conspiracy between the Majorana and Dirac neutrino mass matrices. 
This would appear to be somewhat unnatural in a theory in which Mr and 
N have different origins, as is typically the case in unified models. On the 
other hand, this "Dirac- Majorana conspiracy" might not be unnatural in a 
model in which the same flavor symmetry, and the same small parameter 
characterizing the breaking of that symmetry, controlled the structure of 
both matrices. A good example of such "correlated hierarchies" is the model 

of im 



A very important question is whether 9 a t m can naturally be of order 
unity even if the Dirac matrices are hierarchical and the parameters in Mr 
have no direct relationship to those of N. The answer is yes. In [E3| and 



| 19| interesting examples were found that satisfy these criteria. The specific 
forms given in those papers happen to lead to the SMA solar solution, but 
with some modifications they can also yield a satisfactory LMA solution, as 
we will now see. 

Example 1: The following structure is closely related to that in |2 



N 





( de' ee' 

ge' ae be \m N , M R = [ 1 | m R . 

y he 1 ce 

Here a, h are of order one, e' <C e <C 1, and (e'/e)e _1 C ^4 € e~ l . 
Keeping only the significant terms, the resulting light neutrino mass matrix 
M v = —N T Mft l N is 



/ 0(e' 2 ) O(ee') de'/A \ 2 

M v = - O(ee') a 2 e 2 abe 2 + ee'/A (11) 

V de'/A abe 2 + ee'/A b 2 e 2 + 2fe'/A J mR 

A rotation in the 23 plane by angle 9 = tan -1 (a/6) ~ 1 diagonalizes the 2-3 
block and brings the matrix to the form 



13 



Ml. 



0[e 
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>2\ 



bd 



6- 



If one assumes that e'/A 



e'/A 
e'/A 



ad 

s/a 2 +b 2 
2a(af-be) 

a 2 +b 2 





e'/A 
e'/A 



bd 



Va 2 +b 2 





e'/A\ 



la' 



m 



N 



+ b 2 )e 2 J 



m R 



(12) 



~ e 



! /10, 



then it is apparent that Am^/Am^ ~ 
10 -2 as required for the LMA solution. It is to be observed that the 12 and 
21 elements of this matrix are of the same order as the 22 element. This is 
just what is needed to get the right value of 9 so i for the LMA solution, i.e. 
a value that is of order unity, but not very close to maximal. For example, 
if the 12 element is exactly equal to the 22 element, then tan 2 9 so i = 0.39, 
which is in excellent agreement with the LMA best-fit value given in Q. 
An examination of this matrix reveals that in obtaining the LMA solution a 
crucial role is played by the "cascade hierarchy" form of N. In particular, it 
is important that d be of the same order as e and /, which would not be the 
case if N had a "geometric hierarchy" form. It should also be noted that the 
largeness of the atmospheric angle is also traceable to the cascade hierarchy 
form of N, and specifically to the fact that b is of the same order as a. 



Example 2: The following structure is closely related to that given in [19 



N 



As in the last example, 




m N , 



M, 



R 




m R . 



(13) 



h are of order one, and e' < £ < 1. If one 



also assumes that e 2 /A ^> e' 2 /B 3> 1, then the light neutrino mass matrix 
M v = —N T M^ 1 N takes the form (keeping only the important terms): 



M,j 



0(e' 2 /A) 
gaee'/A + dee' 2 / B 
gbee'/A + dfe' 2 /B 



gaee'/A + dee' 2 / B gbee'/A + dfe' 2 /B 



a 2 e 2 



/A 



e 2 e' 2 



IB 



abe 2 /A + efe' 2 /B 



abe 2 /A + efe' 2 /B 
e 2 /A + f 2 e' 2 /B 



m 



N 



b 2 " 2 



m R 



A rotation in the 23 plane by angle 9 = tan 1 {a/b) 
block and brings the matrix to the form 



1 diagonalizes the 2-3 



Ml, 



0(e' 2 /A) 



d(be-af) n 
Va 2 +b 2 



\ 0(ee'/A) 



IB 
0(e' 2 /B) 



d(be-af) J2 

Va^+U 1 
(be~af) 2 
a 2 +b 2 



e' 2 /B 0(ee'/A) + 0(e' 2 /B) \ 
/B 







(a 2 + b 2 )e 2 /A 



m 



N 



m R 
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The same remarks apply as in the previous example. The largeness of both 
the atmospheric neutrino angle and the solar neutrino angle can be traced to 
the "cascade hierarchy" form of N. Because the 12, 21, and 22 element are 
of the same order, the solar angle is (as required for the LMA solution) of 
order one, but not very close to maximal. The right ratio of mass splittings 
for the LMA solution can be obtained if e /2 /B ~ 10 _1 e 2 /A 

These two examples show that there are reasonable forms or "textures" 
for the mass matrices in the context of the see-saw mechanism that can quite 
naturally yield the LMA solution. However, actual detailed models based on 
these textures have not been constructed. It is also not clear how simple it is 
for the seemingly required "cascade hierarchy" form to arise in the framework 
of grand unified models. Finally, it should be noted that while some forms for 
N and Mr can be identified which would naturally give the LMA solution, 
most of the viable forms give the SMA or VAC solutions, and indeed the 
great majority of see-saw models published in the literature give the latter 
solutions rather than the LMA solution. 

B. See-saw models where 8 atm comes largely from L. 

As we saw in part B. of section 2, there are advantages to models in which 
the large atmospheric angle comes primarily from the charged lepton mass 
matrix L. In particular it becomes easy to reconcile the largeness of 9 atm with 
the smallness of Am 2 ol / Am 2 tm , since they come from different matrices: the 
former from L and the latter from M v . As we also saw, however, having 
a large angle arise from the diagonalization of L raises the question of why 
a large CKM angle does not arise from the diagonalization of the quark 
mass matrices D and U. The answer in the "flavor democracy" models, was 
that the CKM angles are small by a cancellation caused by an approximate 
symmetry. The possibility of a different and very elegant answer arises in 
the context of grand unified models, especially if an SU(5) symmetry plays 
a role in the form of the fermion mass matrices. 

In SU(5) the left-handed (right-handed) charged leptons are in the same 
multiplets as the (CP conjugates) of the right-handed (left-handed) down 
quarks. Since the large mixing angle 6 atm is a mixing of left-handed lep- 
tons, and specifically of left-handed charged leptons in the scenarios we are 
presently considering, it would generally be related by SU(5) to a large mix- 
ing angle for the right-handed down quarks. But such a right-handed mixing 
angle has nothing to do with the observed CKM angles. On the other hand, 
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the small CKM angles are related by SU(5) to small mixings of the right- 
handed leptons, which are irrelevant to neutrino oscillation phenomena. 

These considerations lead naturally to the idea that the matrix L is such 
as to give large left-handed mixings and small right-handed mixings, so that 
Qatm can be large while V c b is small. In other words, L must be highly 
asymmetric or "lopsided" to use the term suggested in [23|. (It should be 
noted that in SU(5), L is related to D T , but not in general to U or N. Thus 
while the lopsidedness of L entails the lopsidedness of D, there is no reason 
to expect N and U to be lopsided, and in the examples we give below they 
are not.) 

Many models have been proposed based on this idea of lopsided mass 

m 



matrices 



These are classified as Type 11(2) in ||. The great majority 
of these models give the SMA solution or the VAC solution to the solar 
neutrino problem. However, it is possible to obtain the LMA solution as well 



26fl. In a lopsided model in which the LMA, LOW, or VAC solution arises, 
the large atmospheric angle can come primarily from the matrix L while the 
large solar angle can come from the matrix M v . This actually has the virtue 
of simplicity, since the form of M v is less constrained than in models where 
it must give rise both to a large 6 a tm and large 9 so i. 

An example of how the LMA solution might arise in a see-saw model with 
lopsided L is provided by the following matrices: 



N 




m R . 



(14) 



/ de' ee' fe' \ 

ge' ae be m N , M 
V he' ce 1 / 

As before, it is assumed that a, h are of order one, and that e' < e < 1. 
Suppose the value of A is such that e 2 , e' <^ ee' /A <C 1. Then, keeping only 
the most important terms, the light neutrino mass matrix has the form: 



( 0{e' 2 /A) 
adee' j A 
\ bdee'/A 



laf 



adee' /A 
2aeee'/A 
I- be)ee'/A + ce 



(af 



bdee'/A 
- be)ee'/A 
1 



ce 



N 



m R 



Here the 23 and 32 elements are small compared to the 33 element, leading 
to a small contribution to 9 atm ; but that is alright, since 9 atm is supposed 
to arise primarily from diagonalizing L in this class of models. As in the 
previous examples, one sees that the 12, 21 and 22 elements are of the same 
order, giving a large, but not nearly maximal, value of 9 so i as required by the 
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LMA solution. To get the right ratio of neutrino mass- squared splittings one 
needs ee'/A ~ 10 -1 . 



4 An 577(5) pattern 

A particularly interesting kind of pattern can arise very simply in the context 
of SU(5) with abelian flavor symmetry. Consider an SU(5) model with a 
U(l) flavor symmetry under which the quark and lepton multiplets have the 
following charge assignments: 10i(2), 10 2 (1), 10 3 (0), 5i(l), 5 2 (0), 5 3 (0), Let 
the breaking of the U(l) flavor symmetry be done by a field \ having U(l) 
charge -1, and an expectation value (%} /Mf [avor = f <C 1. Then the mass 
matrices of the quarks and charged leptons will have the following forms. 



/ 



D 



: 3 e 2 e \ / e 4 e 3 e 2 

e 2 e 1 m D , U ~ e 3 e 2 e | mu, 
e 2 e 1 J \ e 2 e 1 



e 3 e 2 e 2 



(15) 



L ~ | e 2 e e | mo. 

e 1 1 



Note that L and D have the lopsided form. It has been pointed out in several 
papers in the literature that these forms give a very good account of the mass 



ratios and mixing angles of the quarks and leptons [^7 . 

One can see that the quantities m M /m r ~ 1/17, ms/m^ ~ 1/50, and 
V c b — 1/25 all are of order e. Thus, e is roughly of order 1/20. A consistent 
value of e is obtained from the fact that m c /m t ~ 1/400, m u /m c ~ 1/200, 
m e /m^ ~ 1/200, and V u b — 1/300 are all of order e 2 . From the Cabibbo 
mixing and the ratio m d /m s , one would get the somewhat larger value e ~ 
1/5. 

The light neutrino mass matrix M v arises from the see-saw mechanism; 
so to know this matrix exactly it would be necessary to know Mr. However, 
to know merely the order in e of the elements of M v it is not necessary to 
know the U(l) family charges of the right-handed neutrinos at all, since the 
effective mass term in which M v appears involves only the left-handed lepton 
doublets, which are in the 5j. Knowing the U(l) charges of the 5, tells us 
that 
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M v ~ I e 1 1 | m 2 D /m R . (16) 

From the forms of L and M„ it is obvious that the mixing of the 
second and third family neutrinos will get 0(1) contributions from both these 
matrices, thus explaining the largeness of the atmospheric neutrino mixing 
angle. Let us imagine now diagonalizing the 2-3 block of M v to get 



/ e 2 



e m 2 ( ) | m 2 D /m R . (17) 
V c o 

The entry m 2 ( ) would naturally be expected to be of order 1. However, for 
the ratio r = Am 2 oi /Am 2 tm to come out to be of order 10~ 2 , as required by 
the LMA solution, m 2 (o) should rather be of order 1/10. If we accept this 
rather mild fine-tuning, and assume that m 2 (o) ~ 1/10, something interesting 
can be observed, namely that the 12 and 21 elements of M' u are of the same 
order as the 22 element, since e ~ 1/20. Recall that this is just what is 
needed for tan 2 9 so i to come out to be near the best-fit LMA value of about 
0.3 or 0.4. 

This model, then, naturally explains both the value of 9 a tm and the LMA 
value of 9 so i, provided that r = Am 2 oZ /Am 2 im is set to the LMA value. 

Let us now imagine diagonalizing M' v . The rotation needed to eliminate 
the 13 and 31 elements will give a contribution to U e ^ that is of order e, 
quite consistent with the present experimental limit of 0.15. This leaves the 
diagonalization of the 1-2 block. In doing this one may neglect the 11 element, 
since it is of order e 2 . One then finds the simple relations (a) tan2# so z ~ 
2e/m 2 (o), and (b) m 2 = m 2 ( )/(l — tan 2 6 so i)- (Here we have ignored the O(e) 
contribution to 9 so i coming from diagonalizing L, since we are interested in 
large values of 6 so i.) From these relations one can infer roughly what region 
this model gives in the standard log(tan 2 9 so i)-\og( Am 2 ol ) plot. One sees 
from (a) that tan 2 9 so i ~ e 2 (m 2 ( ))~ 2 , and from (b) that Am 2 oZ ~ (m 2 (o)) 2 . In 
other words, in the standard plot the region coresponding to this model lies 
roughly along a line with slope —1 going through the LMA allowed region. 
We shall see shortly, both by a much more careful analytic calculation and 
by a Monte Carlo numerical calculation, that this conclusion is correct. The 
form of (b) tells us something else that is interesting. As the value of the 
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solar angle approaches maximality, i.e. tan 2 6 so i — > 1, the denominator in 
(b) approaches zero. Therefore, to maintain a finite value of Am^ the value 
of 7712(0) must be tuned to be extremely small. Thus, one expects the region 
of greatest probability in this model to fade away as tan 2 9 so i approaches 1. 
This is confirmed by the analytic and Monte Carlo calculations, as we shall 
see. 

We shall now study the predictions of this model in a statistical way, 
much in the spirit of p8|j. Similar statistical analyses have been done in 



several recent papers ||29|| , and our results are consistent insofar as they can 
be compared with theirs. However, our analysis differs in some respects. We 
do not treat e as a free parameter, and seek to find its optimal value for the 
various solar solutions. Rather, we fix e to the value that best reproduces 
the mass ratio m At /m T and then derive the full region of the tan 2 9 so i -Am 2 sol 
plane which results. We also show that by treating the random variables as 
having a Gaussian distribution the statistical predictions of the model can 
be obtained analytically. We also carry out a numerical simulation using 
a non-Gaussian distribution similar to those used in previous analyses and 
show that it agrees remarkably well with the analytic results obtained using 
a Gaussian distribution. 

To carry out the statistical analysis we parametrize the neutrino and 
charged lepton mass matrices as follows: 

/ /e 2 de/V2 ee/V2\ 

M v = de/V2 b c/y/2 m 2 D /m R , (18) 

V ee/V2 c/V2 a J 

and 



L = 0(e 2 ) De Ce m D . (19) 




We will take the unknown order-one parameters a, b, c, d, e, /, A, B, C, 
and D, to be complex random variables whose real and imaginary parts have 
Gaussian distributions with standard deviation a. For example, if a = \a\e l0a , 
then P(a) da = {2ira 2 )~ 1 exp(— |a| 2 /2<r 2 )|a| d|a| dO a . It should be noted that 
we have put factors of l/y/2 in the off-diagonal elements of M v . This is the 
appropriate normalization to use for a symmetric matrix. 

What we want to calculate is the probability distribution P(r, t) dr dt, 
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where r = Am^/Am 2 ^, as before, and t = tan 2 9 so i, given that the order- 
one unknown parameters in the mass matrices have Gaussian distributions as 
described. We will first describe and give the results of an analytic calculation 
of P(r, t), and then present the results of a Monte Carlo numerical calculation 
of P(r,t). 

A very important point in what follows is that if one does unitary changes 
of basis 

(z)- v (z) ° r (I:)^(i)< 

the resulting parameters a' , e', A' , D' , have exactly the same Gaussian 
distributions as the parameters in the original basis. (This would not be 
true without the factors of I/a/2 in M u .) This is one fact that makes the 
analytic calculation tractable using Gaussian distributions. Moreover this 
basis independence is more consistent with the group-theoretical approach 
advocated in p8| . 

The first thing to do is diagonalize L. For our purposes, we need only 
diagonalize the 2-3 block to find m^/m T and the contribution of L to 6 atm . 
This involves multiplying the 2-3 block of L from the right (which in our con- 
vention is a transformation on the left-handed leptons) by a unitary matrix 

ur*(* B ^)(i^i 2 + i5i 2 r 1/2 . 

This eliminates the large element B in L and makes the 33 element become 
\J\A\ 2 + \B\ 2 . The new 22 and 23 elements, which can be written D'e and 
C'e respectively, have the same Gaussian distribution as do A, B, C, and D. 
Consequently, one has that 

(^) ms /K) rms = e\D'\ rms /(^\A\ 2 + \B\ 2 ) rms = e/V2- 

Thus, the most reasonable value to choose for the small parameter, from the 
point of view of lepton physics, is e/ \f2 = m^/m T . 

The first constraint that we shall impose is that the atmospheric neutrino 
mixing comes out to be very close to maximal, as found experimentally. This 
angle gets contributions from the diagonalizations of both L and M v . It 
would seem, then, that we must, in computing P(r,t), take into account 
the random variables in both mass matrices. However, a great simplification 
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occurs because of the use of Gaussian distributions and the resulting basis 
independence of the probability distributions. A little thought shows that 
one can compute P(r, t) in a basis where the contribution to the atmospheric 
neutrino mixing coming from L has some fixed value, and the result will not 
depend on that value. Thus the parameters in L are irrelevant to P(r, t). It is 
simplest in practice to choose the basis where the entire atmospheric neutrino 
mixing comes from L. A further simplification is achieved by neglecting the 
parameters e and /. The parameter e comes into calculating U e s (which is 
predicted to be of order e, and so consistent with the experimental bound 
| C/e3 1 < 0.15), but has a negligible effect on 9 so i, 9 a tm, and the neutrino masses. 
The parameter / is multiplied by e 2 , and so is negligible also. Finally, one 
can choose a basis where the parameters a, b, and d are real. That means 
that the only parameters that come into the calculation of P(r,t) are \a\, \b\, 
\c\, \d\, and 6 C = argc. From now on, we shall drop the absolute value signs 
and denote \a\, for example, simply by a. 
One begins, then, with a matrix 

/ de/V2 \ 

M v = de/V2 b ce idc /V2 m 2 D /m R , (20) 
V a ) 

and a probability distribution 

P(a,b,c,9 c ,d) = ^e-^ 2 ^ 2 ^. (21) 

The first step is to diagonalize the 1-2 block of the matrix given in Eq. 
(20), which gives tan2# so ; = s = \f2dejb. This allows the elimination of the 
random variable d in favor of the measureable parameter s or equivalently 
t. The 11 and 22 elements of the matrix then become |(1 — \/l + s 2 )b and 
1(1 + + s 2 )b respectively. The latter quantity we shall denote as b' . The 
next step is to impose the atmospheric angle constraint. Since we are working 
in a basis where the contribution to this angle from M v is vanishingly small, 
the imposition of this constraint sets the parameter c to zero. More precisely, 
if one requires that the (complex) contribution to the atmospheric angle 
from M v have magnitude bounded by some arbitrary small cutoff A aim 1, 
the condition on c becomes (c/y/2)/(a — V) < A atm . This means that the 
integration over dcd# c in the probability distribution can be done, yielding 
/cdcd# c = n(y/2(a - b')A atm ) 2 . The only remaining variables are then a, 6, 
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and s. The random variable a can be eliminated in favor of the measurable 
ratio r of mass-squared splittings using the relation r = (b 2 \/l + s 2 ) / (a 2 — b' 2 ). 
It is a very good approximation here to replace b' by b, since for the whole 
region of interest either s or r is very small as we shall see. After all these 
steps, one is left with a probability distribution P(r,b,s). The final step 
is simply to integrate over the random variable b. Since this integral is a 
Gaussian it is easily done. The final result is 



P(r, s) drds = N 



rs dr ds 




VlT 



1 



2r 



+ 



rs 



-I 4 ' 



(22) 



Vl + s 2 ' 2eVl + s 2 _ 

where N is just a normalization constant. Changing variable from s = 
tan26L n / to t = tan 2 6L„ ; one finds that 



P(r,t) drdt = N 



2r dr dt 




2r 



t 



2rt 



(23) 



One can see the qualitative behavior of this function rather easily. The 
crucial term is the one containing e 2 in the denominator. This term forces the 
product rt to be of order e 2 . This is consistent with what we argued above, 
namely that the region of greatest probability in this model has rt ~ constant, 
i.e. a line of slope -1 in the log(tan 2 6' so ;)-log(Am 2 oi ) plane. Moreover, we 
see that as t — ► 1, the product rt is forced to be less than or of order 
e 2 (l — t 2 ) — > 0, so that the probability is suppressed for t = 1. 

In Figure [I], we give a contour plot of the probability function just com- 
puted analytically, and compare it to the results of a Monte Carlo calculation. 
For the Monte Carlo calculation we used the forms in Eqs. (18) and (19), with 
e/\^2 = m^/rrir, but assumed that the magnitudes of the complex random 
variables, instead of having a Gaussian distribution, had constant probabil- 
ity in the interval 0.5 to 2.0 and zero probability outside that interval. The 
phases of the complex variables were also treated as random numbers and 
were varied from to 2ir. We then diagonalized randomly generated matri- 
ces to obtain the corresponding MNS mixing matrices Umns and neutrino 
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masses and analyzed the results by imposing the conditions sin 2 2# otm > 0.9 
and \U e3 \ < 0.15. These conditions reduced our initial set of 50,000 data 
points to 20,860 that were compatible with both the CHOOZ and the atmo- 
spheric neutrino experiments. The points that passed the cuts are given in 
Fig. |. 

One can see from the excellent agreement between the analytic and Monte 
Carlo results evident in Fig. [I] that the exact form used for the probability 
distributions of the random variables themselves makes little difference. This 



has also been found in other papers [gq, 29 . One point that should be noted 



is that since Figure |T] is a log-log plot, the correct thing to plot and what has 
been plotted, is P(logr, log£) ~ P(r,t)rt. 

In Figure we have taken the slice r = 1.4 x 10~ 2 , which comes from using 
the best-fit values from experiment, and plotted the resulting P(logi) against 
a normalized, binned distribution. The binned distribution has been obtained 
by counting the number of data points in the strip logr = log(1.4x 10 _2 )±0.1, 
the width of one bin being 0.2, and the normalization has been carried out 
with respect to the maximum of P(logt). Note that the most probable value 
for tan 2 6 so i is about 0.1, with a very substantial part of the area under the 
curve being in the region [0.2, 0.8] preferred by the LMA solution global fits. 

The one weakness of this model is that it does not explain why r = 
^ m L;/^ m atm ~ 1-4 x 10~ 2 . From Figure |I]one sees that a value of 10 _1 for 
this ratio is near the peak of the probability distribution P(r, t). However, the 
same figure shows that a value of 1.4 x 10 -2 is near the edge of the preferred 
region, and so requires a mild fine-tuning. However, once r is constrained to 
be the right value, the value of tan 2 9 so i needed for the LMA solution emerges 
quite naturally, as can be seen from Fig. |^. The atmospheric mixing angle 
is, of course, also naturally explained. In Figure [l| the best fit values of the 
LMA, SMA, and LOW solutions are indicated by dots. One sees that this 
SU(5) lopsided model (with e/y/2 = m^/m T ) naturally prefers the LMA 
solution over the others. 



5 Conclusions 

One can see from the foregoing that it is considerably easier to build satisfac- 
tory models of the VAC, LOW, or SMA type than of the LMA type. That is 
reflected in the models that have actually been constructed in the literature. 
One problem is that in many models which predict large solar mixing angle, 
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notably inverted hierarchy schemes and flavor democracy schemes, this an- 
gle tends to come out very close to maximal. They do not naturally explain 
why tan 2 9 so i should come out in the range 0.2 to 0.8 preferred by the data. 
Other non-see-saw schemes, such as SUSY with R-parity breaking and single- 
right-handed-neutrino-dominance (SRHND) models, tend to predict a value 
of Am 2 so J Am 2 atm significantly less than that preferred by the LMA solution. 
While the LMA value of this ratio can be fit, it is not really explained. 

The situation seems more promising for the see-saw approaches, although 
here also the great majority of published models give the small angle or 
vacuum solar solutions. We showed that certain fairly simple textures exist 
that would naturally reproduce the neutrino masses and mixings required by 
the LMA solution. Whether these textures can be implemented in simple 
models remains to be seen. 

One of the few existing schemes that shows a natural preference for the 
LMA solution is the lopsided SU(5>) model studied in Section 4. The value 
of Am 2 sol / Am 2 atm requires a mild fine-tuning, but given that, both the atmo- 
spheric angle and the LMA value of the solar angle emerge quite naturally. 
We studied the predictions of this model in a statistical way, and found that 
by using Gaussian distributions the analysis could be carried out very simply 
and accurately by purely analytic means. We believe that the same methods 
should be applicable to many other models. The advantage of such statisti- 
cal analyses is that they allow one to estimate in a somewhat objective and 
quantitative way how "fine tuned" models must be to reproduce the data. 
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Figure 1: Contour plot of the normalized probability distribution P in 
log(tan 2 9 so i)-\og(Aml ol / Aml tm ) plane with the contour values 0.002, 0.02, 
0.06, 0.1, and 0.14 superimposed on the numerically generated distribution 
of points. Large dots represent best-fit values for LMA, SMA, and LOW 
solar neutrino solutions, as indicated. 
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Figure 2: Normalized probability distribution P (solid line) as a function 
of log(tan 2 9 so i) for best-fit LMA solution value Am 2 sol /Aml tm = 1.4 x 10 -2 
plotted against the normalized, binned distribution (dashed line). 
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